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Abstract 

We calculate flavor-changing neutral-current (FCNC) processes involving top-quark 
production via anomalous tqV couplings at the Tevatron and HERA colliders. We cover 
the FCNC processes pp — > tZ, pp — ► tj, pp — ► tt, and ep — ► el. We go beyond leading order 
and include the effects of soft-gluon corrections through next-to-next-to- leading order. We 
demonstrate the stabilisation of the cross section with respect to the variation of QCD 
scale, and we investigate the reach of the Tevatron and HERA colliders. 



1 Introduction 



The top quark, the heaviest known fermion, occupies a unique place in the search for new 
physics beyond the Standard Model (SM). This is the only fermion with a mass of about the 
scale of the electroweak symmetry breaking (EWSB). Therefore the study of its properties 
and their possible deviations from SM predictions could shed a light on the mechanism of 
EWSB. Top quark physics can probe various physics beyond the SM: anomalous gluon-top 
quark couplings pQ, anomalous Wtb couplings [2], new strong dynamics 3j, flavor- changing 
neutral-currents (FCNC) [I], R-parity violating SUSY effects [5], CP- violation effects |3, and 
effects of Kaluza-Klein excited P^-bosons [7j. 

In particular, the study of FCNC couplings involving the top quark is well motivated. The 
effective Lagrangian involving such couplings of a t, q pair to massless bosons is the following: 

AC eff = jK tqV eia^qF$ u + h.c, (1.1) 

where n tq v is the anomalous FCNC coupling, with q a u- or c-quark and V a photon or Z-boson; 
Fy are the usual photon/Z-boson field tensors; = {i/2){^^y u — 7^7^) with 7^ the Dirac 
matrices; e is the electron charge; and A is an effective scale which we will take throughout this 
paper to be equal to the top quark mass, which we denote by m. In the SM these operators can 
be induced through higher-order loops, however these effects are too small to be observable [8J. 

In some models involving new physics FCNC top-quark couplings could appear at tree- 
level and, therefore, lead to large FCNC effects. Such enchancements of FCNC couplings 
involving top quarks could appear in two-Higgs doublet models and supersymmetric models [SJ 
as well as in multiple-Higgs doublet models [TU]. Models with FCNC coupled singlet quarks, 
dynamical electroweak symmetry breaking, or compositeness, can further enchance FCNC top- 
quark couplings 11 . 

The presently operating colliders such as the Tevatron and HERA have a nice opportunity 
to probe FCNC interactions in the top-quark sector. In order to establish accurate limits on 
FCNC couplings, experiments need accurate predictions of cross sections for FCNC processes. 
We have shown in that at HERA uncertainties for the tree-level FCNC cross section ep — > et 
can be very big, as the cross section can vary by a factor of two due to the choice of the QCD 
scale. This fact unavoidably requires involving the higher-order corrections for the stabilization 
of the FCNC cross section. 

In this paper we study some principal FCNC processes both for the Tevatron and HERA 
colliders involving t — u — V couplings with V a photon or Z-boson ^3] • The effects of 7 — Z 
interference are also taken into account in our studies. We include soft-gluon corrections at next- 
to-leading order (NLO) and next-to- next-to-leading order (NNLO). We present the stabilization 
of the NNLO cross sections with respect to the variation of the renormalization/factorization 
scales. 

In Sections 2 and 3 we study FCNC processes of the associated top-quark production with 
Z-boson and photon at the Tevatron, respectively. Section 4 is devoted to the uu — > tt process 
of same-sign top quark production at the Tevatron. Section 5 extends our previous study on 
single top-quark production at HERA, and we now consider this process at NNLO level as well 
as study both t — u — 7 and t — u — Z vertexes and take into account 7 — Z interference. Finally 
in Section 6 we draw our conclusions. 

We note that contributions involving the charm quark via t — c — V couplings are greatly 
suppressed because of the small charm parton densities. In each section we make a note about 
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the contribution from charm quarks. We also note that we don't consider ti production and tq 
or tq production in this paper because the Standard Model cross sections for these processes 
overwhelm any FCNC contributions. 

In the present paper we use the notation a for parton-level differential cross sections, while 
a denotes the hadronic cross sections. We also define a kinematical variable S4, for each 
subprocess, that goes to zero at threshold. The effects of the soft-gluon radiation appear 
in the form of logarithmic "plus" distributions with respect to S4 of the type 



r hr fe (s4 / /m 2 



s 4 



k < 2n - 1 



;i-2) 



at n-th order in the QCD coupling a s . These plus distributions are the remnants of cancellations 
of infrared divergences between soft and virtual contributions to the cross section. We define 
the leading logarithms (LL) with k = In — 1 and the next-to-leading logarithms (NLL) with 
k = 2n - 2. 

The plus distributions are defined by their integral with any smooth function /, such as 
parton distributions, as 
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They provide significant and dominant contributions to the cross section near threshold, where 
there is limited phase space for real gluon emission, as has been shown for many Standard Model 
processes, including top production [Hj, direct photon production [TH], jet production [IE] , and 
P^-boson production J7j. A unified approach for the calculation of these soft-gluon corrections 
for various processes in hadron-hadron and lepton-hadron colliders has recently been presented 
in Ref. [H]. We follow that reference in calculating soft-gluon corrections to FCNC processes 
at NLO and NNLO with NLL accuracy, i.e. keeping LL and NLL. We work in the MS scheme 
throughout the paper. 



2 gu — > tZ 

We begin with FCNC tZ production at the Tevatron. For the process g{p g ) + u{p u ) — > t(p t ) + 
Z(pz), we define the kinematical invariants s = {p g + p u ) 2 , t = (p g — p t ) 2 , u = {p u — p t ) 2 , and 
S4 = s + t + u — m 2 — m|, where mz is the Z-boson mass and m is the top quark mass. Note 
that near threshold, i.e. when we have just enough partonic energy to produce the tZ final 
state, S4 — > 0. The respective tree-level Feynman diagrams are shown in Fig. ^ 




Figure 1: Tree- level Feynman diagrams for the process gu — > tZ. 
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The differential Born cross section is 



d 2 a B u ^ tz 
dt du 



rpgu^tZ 
r B 



S{s 4 ) 



(2.1) 



where 
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where a = e 2 /(4ir) and for brevity we define Kz = Ktuz- 
The NLO soft-gluon corrections for gu — > tZ are 
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Here cf ^ tz = 2(C F + C A ), where C F = (N 2 - 1)/(2N C ) and C A = N c with N c = 3 the 
number of colors, and 
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where fj, F is the factorization scale, and we have defined r £^ u ^' tz as the scale-independent part 
of C2 U_> * Z . The term Rer'^ denotes the real part of the one- loop soft anomalous dimension, 
which describes noncollinear soft-gluon emission jTHj. A one-loop calculation gives 
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c 



gu-^rtZ 



C F ln 



—t + mi 



C A \n 



—u + m z 



3^, /?o 
— o p — — 
4 4 



/4 



/?0 



In ^ +-ln ^ 



Ma' 



(2.6) 

where lir is the renormalization scale and (3q = (HCa — 2ny)/3 is the lowest-order /3 function, 
with rif = 5 the number of light quark flavors. Note that cf u ~^ tz represents the scale-dependent 
part of the S(s4) corrections. We do not calculate the full virtual corrections here. Our calcu- 
lation of the soft-gluon corrections includes the leading and next-to-leading logarithms (NLL) 
and is thus a NLO-NLL calculation. 
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We next calculate the NNLO soft-gluon corrections for gu — > tZ: 

d 2 af a ^ tz _ ^ tz a^(jj^) J 1 / gu ^ tz \2 [ ln 3 (s 4 /m 2 ) " 
dtdu vr 2 i 2 V c 3 
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where (2 — tt 2 /6 and £3 = 1.2020569.... We note that only the leading and next-to-leading 
logarithms are complete. Hence this is a NNLO-NLL calculation. Consistent with a NLL 
calculation we have also kept all logarithms of the factorization and renormalization scales in 
the [ln(s 4 /m 2 )/s 4 ] + terms, and squares of scale logarithms in the [l/s 4 ] + terms, as well as ( 2 
and (3 terms that arise in the calculation of the soft corrections. For relevant details in the 
related process of Standard Model top production see Ref. [T4"] . 

We now convolute the partonic cross sections with parton distribution functions (PDF) 
to obtain the hadronic cross section. For the FCNC hadronic cross section p(p a ) + p{Pb) — * 
t(p t ) + Z{pz) we define S = (p a + Pb) 2 , T = (p a — p t ) 2 , and U = (p& — p t ) 2 , and note that 
Pg — %aPa, Pu — %bPbi where x denotes the momentum fraction of the hadron carried by the 
parton. The hadronic cross section is then given by 
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(2.9) 



x, 



\s 4 



m 



2 y *' 2 

m| — Xb(U — m 2 )] I (xf,S + T — m 2 ), and <p(x) are the parton distributions. 
In our calculations we are using the MRST2002 NNLO parton distribution functions j20] 
with the respective three-loop evaluation of a s . For all three — Born, NLO-NLL, and NNLO- 
NLL — we use the same MRST2002 NNLO PDF. We also use \ip = /ir for our numerical 
results. For reference point we choose n z =0.1 and fi = \ip = fiR = m = 175 GeV. In our paper 
we quote the theoretical error related to the conventional variation of the scale between m/2 
and 2m. 

We present our numerical results in Fig. El The left frame presents the cross section for 
the process gu — > tZ versus top-quark mass at the Tevatron with y/S = 1.96 TeV, at Born, 
NLO-NLL, and NNLO-NLL orders. Note that the first-order corrections, A<tnlo-nll, and the 
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Figure 2: The Born, NLO-NLL, and NNLO-NLL cross sections for gu — > tZ in pp collisions 
with y/S = 1.96 TeV and k z = 0.1. Left: cross section versus top-quark mass at scale /i = m; 
right: cross section versus the scale for m = 175 GeV. 

second-order corrections, Aa NNLO _ NLL , are positive for /i = m. There is clear stabilisation of 
the cross section due to the Ao- NLO -nll and Act nnlo—nll corrections over a large range of 
n/m presented in the right frame of Fig. |21 At the reference point, fi = m = 175 GeV and 
Kz = 0.1, we have the following results: 



a 



B^n = 61 ft>, ^NLO-NLL = 20 ^NNLO-NLL = 6 fb. 



Taking into account the theoretical error due to the scale variation between m/2 and 2m one 
thus has for the total NNLO-NLL cross section: 

gu^tZ _ 07+2 fi 

a NNLO-NLL — °'-3 iu - 

Since the total cross section is proportional to k\ it is straightforward to rescale it for any 
given value of the anomalous FCNC coupling. 

For completeness we would like also to quote the cross section for the case of n tc z coupling. 
Since the initial c-quark is coming from the sea, and thus the charm density is quite small, the 
total cross section for gc — > tZ at the Tevatron, with [l = m — 175 GeV and assuming the 



anomalous coupling is again 0.1, is about 40 times smaller: <7 9 ^nlo 



NNLO-NLL 



2.4 fb. 



One should also notice that we present all cross sections for top-quark production. For 
experimental studies one should also include the contribution from anti-top-qu&rk production 
which will double the total cross section for all processes at the Tevatron we study here. 



3 gu — > £7 

We continue with FCNC £7 production at the Tevatron. For the process g(p g ) + u(p u ) 



t{p t ) +7(^7), we define the kinematical invariants s = (p g +p u ) 2 , t = (p g —pt)' 



u 
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and s 4 = s + t + u — m 2 , where again m is the top-quark mass. Note that near threshold s 4 — > 0. 
The respective tree-level Feynman diagrams are shown in Fig. El 
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Figure 3: Tree-level Feynman diagrams for the process gu —>■ £7. 
The differential Born cross section is 
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The NLO-NLL corrections for gu — > £7 are 
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with r '5 defined in Eq. (12. 5 J) , and 
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We note that the Cj coefficients are similar to those for the — > £Z process. This is because 
the QCD content of the two processes is the same. 
The NNLO-NLL corrections for gu — > £7 are 
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Again, they are of the same form as for the gu — > £Z process. 
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Figure 4: The Born, NLO-NLL, and NNLO-NLL cross sections for gu — > £7 in pp collisions 
with \/£ = 1.96 TeV and k 7 = 0.1. Left: cross section versus top-quark mass at scale \x = m; 
right: cross section versus the scale for m = 175 GeV. 



For the FCNC hadronic cross section p{p a ) + p(Pb 
T = (p a — p t ) 2 , and U = (pb — Pt) 2 , and note that j 
section is then given by 

rm 2 +m 2 S/(T-m 2 ) 



t(pt)+l(p 7 ) we define S = (p a +p b ) 2 , 
■ x a Pa, Pu = XbPb- The hadronic cross 
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with x a = [s 4 — m 2 — x b {U — m 2 )] / (x b S + T — m 2 ). 

Our numerical results are presented in Fig.HJ The left frame presents the cross section for the 
process gu — > £7 versus top-quark mass at the Tevatron, with v^S* = 1.96 TeV, at Born, NLO- 
NLL, and NNLO-NLL orders. One can see that the NLO-NLL correction is positive at \x = m 
but the NNLO-NLL correction is negative. The right frame presents the scale dependence of 
the cross section which again is stabilized when corrections are added. At /i = m = 175 GeV, 
the Born and NNLO-NLL results happen to be very close to each other because of the negative 
NNLO-NLL correction. For our reference point, \i — m— 175 GeV, « 7 = 0.1, we have: 
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This gives a NNLO-NLL total cross section 



& NNLO-NLL — yd -ll iu 

with theoretical error due to the scale variation between m/2 and 2m. 

Since the total cross section is proportional to k 2 it is straightforward to rescale it for any 
given value of the anomalous FCNC coupling. 

For completeness we would like also to quote the cross section for the case of K tcy coupling. 
Since the initial c-quark is coming from the sea, the total cross section at the Tevatron for the 
process gc — > tj, at n — m — 175 GeV and assuming the coupling is again 0.1, is about 35 
times lower then the one for gu — > tj: ct^nlo-nll = 

2.7 fb. 

Again we note that we present all cross sections for top-quark production. For experimental 
studies one should also include the contribution from anti-top-quark production which will 
double the total cross section at the Tevatron. 



4 uu — > tt 

We now turn to FCNC same-sign top-quark production at the Tevatron. For the process, 
u (Pu a ) + u i.Pu b ) —* t(Pi) + tfa), we define the kinematical invariants s = (p Ua + p Ub ) 2 , t = 
{Pu a ~ Pi) 2 ; t\ — t — m 2 , u = (p Ub — pi) 2 , u\ = u — m 2 , and = s + t\ + «i, with m the 
top-quark mass. At threshold — > 0. The respective tree-level Feynman diagrams are shown 
in Fig. El 
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Figure 5: Tree-level Feynman diagrams for the process uu — > tt. 



The color structure of the hard scattering is somewhat complex since we have two top quarks 
in the final state and two up quarks in the initial state. This is in contrast to all the other 
processes discussed in this paper which have simple color structure. We thus decompose the 
tree-level amplitudes in terms of color tensors using the methods and techniques detailed in 
Ref. JH]- As we will see the soft anomalous dimension is now a matrix in color space. 

For the process at hand we choose the color basis consisting of singlet exchange in the t 
and u channels, c\ = 5 a i^62 and c 2 = 5 a2 5&i, where a (b) is the color index for the w-quark with 
momentum p Ua (p Ub ) and 1 (2) is the color index for the top quark with momentum p\ (p 2 ). Of 
course the physical cross section is independent of the specific choice of color basis. We write 
the Born cross section as a trace of the product of a "hard" function, H, that describes the 
short-distance hard-scattering, and a "soft" function, S, that describes noncollinear soft gluon 
emission |14 [ ITT ^ ITH]. Both H and S are 2x2 matrices in the color basis ci, c 2 . 

The hard matrix at lowest order, calculated by projecting the tree-level amplitude onto the 
color basis, is given by 
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The soft matrix at lowest order is given by 
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Note that the matrix elements of are simply given by S^f = Tr[c*Cj]. 



Then the differential Born cross section is 
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The NLO-NLL corrections for uu — > tt are 
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is again the scale-dependent part of the £(54) terms. Also 
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is the part of the NLO-NLL corrections not proportional to the Born cross section; this is 
again because of the complex color structure of the hard scattering. Here V g is a 2 x 2 soft 
anomalous dimension matrix, calculated at one loop, with elements 
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where Lg = [(1 - 2m 2 /s)//3} [ln((l - /3)/(l + #)) + rri], with 0=^1- Am 2 /s. Note that all 
imaginary parts cancel out in Eq. (J4.12j) . 
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Figure 6: The Born, NLO-NLL, and NNLO-NLL cross sections for uu — > tt in pp collisions 
with y/S = 1.96 TeV and k 7 = Kz = 0.1. Left: cross section versus top-quark mass at scale 
fj, = m; right: cross section versus the scale for m = 175 GeV. 



The NNLO-NLL corrections for uu — > tt are 



T-IUU- 
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-4C F ln(4 



777/ 
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(4.14) 

For the FCNC hadronic cross section p(p a ) + p(p&) — ¥ t(Pi) + t(p2) we define S = (p a + Pb) 2 , 
T = (p a — px) 2 , and U = (pb ~ Pi) 2 , and note that p Ua = x a p a , p Ub = X},Pb- The hadronic cross 
section is then given by 



a 



FCNC 



VP- 



dT 



(m 2 -T)/(S+U-m 2 ) 



-(S/2)(l-y/l-4m*/s) 
-(5/2)(l +v /l-4mVs) 

fX b (S+U-m 2 )+T-m 2 

dxb 



m 2 +m 2 S/(T-rn 2 ) 



S-T+2m 2 



dU 



dsi 



X a Xb 



XbS + T 



<p(x a ) 4>{x b ) 



d 2 a v 



dt du 



(4.15) 



with x a = [s 4 — Xb(U — m 2 )} / (x^S + T — m 2 ). 

Our numerical results are presented in Fig. El The left frame presents the cross section for 
the process uu — > tt versus top-quark mass at the Tevatron, with y/S = 1.96 TeV, at Born, 
NLO-NLL, and NNLO-NLL orders. One can see that the /S.gmlo-nll and Actnnlo-nll 
corrections are always positive for this case for \x = m. The stabilisation of the NNLO-NLL 
cross section with respect to scale variation is even more pronounced for this process compared 
to the two previous cases as one can see in the right frame of Fig. El At \i — m — 175 GeV, we 
have for = Kz — 0.1: 



a 



At 



Born 



1.36 fb, Act 



«i->tt 

NLO-NLL 



0.28 fb, 



A UU—*tt 

LAO NNLO-NLL 



0.10 fb. 



Taking into account the theoretical error due to the scale variation between m/2 and 2m one 
has 

hO.OO 



uu—*tt 
^ NNLO-NLL 



1 74+U.UU n 

L - ' ^-0.02 iu - 



We can investigate the magnitude of the cross section if one of the anomalous couplings is 
set to zero. If we use « 7 = 0.1 and set Kz = the NNLO-NLL cross section is 0.50 fb, while if 
we set k 7 = and use Kz = 0.1 the NNLO-NLL cross section is 0.38 fb. 

We also note that we can have contributions from anomalous couplings involving the charm 
quark. The total cross section for the process cc — > tt at the Tevatron, using n tcl = n tcZ = 0.1, 



at n 



m 



175 GeV, is: a 



*tt 

NNLO-NLL 



0.0044 fb. It is almost 3 orders of magnitude lower 



then the cross section for the uu — ► tt process since both c-quarks originate from the sea. 

The total cross section at the Tevatron for the mixed process uc — > tt, using the value 0.1 
for all anomalous couplings and /i = m = 175 GeV, is: &Wnlo-nll = 0-62 fb. 



12 




Figure 7: Tevatron reach in kz — ft 7 plane for the process uu — > tt for cross section levels 15, 
30, and 45 fb. 



Also, we note that tt FCNC production at the Tevatron has an equal cross section to tt 
production, so if it is included the total cross section is doubled. 

Processes like same-sign top-quark production play a special role since tt final states give 
rise to same-sign leptons (SSL). We have estimated that potential uu — > bbW + W+ background 
can be safely neglected. Assuming 10 fb" 1 integrated luminosity and 1% efficiency (including 
b — tagging and top-quark decay branchings) one should require 30 fb of signal in order to 
observe 3 signal events (to exclude signal at 95% CL according to Poisson statistics under the 
assumption that the background is negligible). We use this cross section as a rough idea on the 
order of the signal cross section to which Tevatron will be sensitive. Since the process uu — ► tt 
depends on two parameters, we present 15-30-45 fb contour levels of the uu — > tt cross section 
in k 7 — Kz plane in Fig. This figure, which gives a rough idea about the Tevatron reach, 
shows that the Tevatron sensitivity is just slightly better for the k 7 coupling. 

At this point, it is worth mentioning the present experimental limits on the k 7 , kz couplings. 
First, a limit on FCNC top-quark couplings has been established by the CDF collaboration in 
terms of the limit on the branching ratio of FCNC top-quark decay [21]: Br(t —>■ q'j) < 3.2%, 
Br(t — > qZ) < 33% which corresponds to k 1 < 0.14, Kz < 0.69 limit in the notation of 
Eq. (jl.lj) . A more stringent limit on k 1 < 0.10 has been established by the L3 collaboration [2*2"] . 
Finally, the best limits on k 1 have been established by HI (k 7 < 0.090) [2HJ E] and by ZEUS 
(k 7 < 0.058) (2Sll2Sll2Zj collaborations. All the limits quoted above correspond to the notation 
of Eq. (|l.lj) . One can see that in light of the present experimental limits, the uu — > tt process 
is unlikely to be observed at the Tevatron unless uu —>■ tt production involves chromomagnetic 
flavor-changing current, which is beyond the scope of the present paper. In order to estimate 
Tevatron sensitivity to pp — > tZ, pp —>■ £7 processes one should study in details all potential 
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backgrounds which is, again, beyond the scope of our paper. 



5 eu — > et 

The last process we present is single-top production at HERA |23 123 123 123 E| We have 
previously studied this process in Ref. [T2j where we calculated the NLO soft corrections with 
tq'j couplings. Here we extend our previous calculation by also including tqZ couplings and 
NNLO-NLL soft corrections. 

For the process e(p e ) + u{p u ) — ► e(pf) + t(p t ), we define the kinematical invariants s = 
(Pe +Pu) 2 , t = (p t — Pu) 2 , u—{pt — p e ) 2 , and = s + t + u — m 2 — 2m 2 , with m the top-quark 
mass and m e the electron mass. Note that near threshold S4 — > 0. The respective tree-level 
Feynman diagrams are shown in Fig. |HJ 

e ► t — > e e ► -j — ► e 

u ► <— » — t u ► <— ► — f 



Figure 8: Tree-level Feynman diagrams for the process eu — *> et. 
The differential Born cross section is 



j2 ~eu—>et 
a a B 

dt du 



T7<eu — >et 
B 



(5.1) 



where 



7B 

eu^ret 



TC a 



m 2 (s — m 2 ) 2 



{-t [2m 4 + m 4 - 2s 2 + (2s + t)(2s - m 2 - 2m 2 )] - 2m 2 m 4 } 



5 ft' 



4:KryKz(l — 4 sin 2 6\ 



.,, (. / ! . i:<]i. "\y ) K z(l ~~ ^ sm2 @w + 8 sin 4 6\y) 



t 2 sin 6 W cos 6*jy £ (t — m|) sin 2 0^ cos 2 0yp (i — m|) 2 



K 2 7 m^-2m 4 + tm 2 + t 2 



sin W cos2 #w (t — rn 2 z ) 2 



(5.2) 



with 9w the weak-mixing angle, k 7 = K tu <y, and = ftt«.z. 
The NLO-NLL corrections for eu — > et are 



d 2 a {1) 



eu^et 



dt du 



w eu^et a s{.^ 2 R ) 



_eu— *et 



7T 



ln(s4/m 2 
s 4 
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1 L-2 



1 

S4 



5(s 4 ) 



(5.3) 



Here c^^ 6 * 



2C F , 



„eu—>et 



2Tf - C F 



+ In 



\ m 2 



(5.4) 
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where = Cp ln[(m 2 — t)/(y/sm)] is the one-loop soft anomalous dimension and Tr 
the scale-independent part of c| M ^ e ', and 



e,u— *et 



3 , 
-- + ln 

4 



—u + mi 



is the scale-dependent part of the <5(s 4 ) terms. 
The NNLO-NLL corrections for eu —>■ et are 



CVln 
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2 lb 



(5.5) 
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(5.6) 



For the FCNC hadronic cross section e{p e ) + p(p P ) —> e(Pf) +t(pt) we define S = (p e +p P ) 2 , 
T — (pt — Pp) 2 , and U = (p t — p e ) 2 , and note that p u = xp p . The hadronic cross section is then 
given by 



a 



FCNC 
ep^et 



(S) 



Tmax 



dT 



x 



S + T 



S-T+m 2 +2m 2 e 
(j)(x) 



dU 



S+T+U-m 2 -2m 2 



dSi 



d 2 a, 



eu—ret 



77T 



dt du 



(5.7) 



with 



rji _ 2 

X max — TTh — 



(g - K) 

2S 



S + m — m e =f y(S - m 2 — m 2 ) — Am 2 m 2 



(5.f 



U max = 2m 2 + m 2 — T — S{m 2 - T)/(S - m 2 ) - (S - m 2 e )m 2 /{m 2 — T) , and x = (s 4 + m 2 - 
t/)/(S + T-m 2 -m 2 ). 

Our numerical results are presented in Fig. El The left frame presents the cross section for 
the process eu — > et versus top-quark mass at HERA, with v^S* = 318 GeV, at Born, NLO- 
NLL, and NNLO-NLL orders. One can see that the Act nlo—nll corrections are positive while 
the Act nnlo—nll corrections are negative for this case for fj, — m. The stabilisation of the 
NNLO-NLL cross section with respect to scale variation is shown in the right frame of Fig. El 
At jj, — m — 175 GeV, we have for k 7 = Kz = 0.1: 



a 



eu^et 
Born 



0.56 pb, Aa™Eo-NLL = O.H pb, Aa 



'nnlo-nll — 0-03 pb. 



15 



eu -> et 



eu -> et 




m (GeV) n / m 

Figure 9: The Born, NLO-NLL, and NNLO-NLL cross sections for eu — > et at HERA with 
Vo = 318 GeV and k 7 = k z = 0.1. Left: cross section versus top-quark mass at scale /x = m; 
right: cross section versus the scale for m = 175 GeV. 



Taking into account the theoretical error due to the scale variation between m/2 and 2m one 
has 

" NNLO-NLL — u - D4 -0.04 P D - 

We note that almost all of the contribution comes from the k 7 coupling. If we use « 7 = 0.1 
and set Kz = the cross section remains practically unchanged (i.e. 0.64 pb at NNLO-NLL), 
while if we set k 7 = and use Kz — 0.1 the NNLO-NLL cross section is only 0.0018 pb. 

In Fig. El we present contour levels for the eu — > et process in k 7 — kz plane. One of 
the contours is at 0.55 pb level which is the present HERA limit reported by HI [23 121] 
collaboration. The second contour is for 0.23 pb level which corresponds to the limit reported 
by ZEUS |2Z] collaboration. The third contour is for 0.1 pb level anticipating the improvement 
of the experimental limits for higher luminosities of the coming HERA II stage. Fig. clearly 
demonstrates that the HERA collider (in contrast to the uu — ► tt process at the Tevatron) is 
much more (about a factor of 40) sensitive to the k 7 coupling than to the Kz one. 

We also note that we can have contributions from anomalous couplings involving the charm 
quark. Using K tcJ = K tcZ = 0.1 the cross section at HERA for ec — > et, with \i = m = 175 GeV, 
is Vnnlo-nll = 0.0020 pb, which is low due to the small sea charm-quark densities. 

In the case of et production, involving the anti-top, the cross section is small, o'n'nlo-nll = 
0.0079 pb at /x = m = 175 GeV, and thus asymmetrical to et production, in contrast to the 
other processes which we considered at the Tevatron. This is because at HERA the valence 
up-quark density can not contribute to et production. 
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Figure 10: HERA reach in Kz — k 7 plane for the process eu —>■ et for cross section levels 0.55, 
0.23, and 0.1 pb. 



6 Conclusions 

We have studied several FCNC processes at the Tevatron and HERA colliders involving top 
quark production via anomalous tqV couplings, with q an up or charm quark and V a photon or 
Z-boson. The partonic processes involved are gu — > tZ, gu — > tj, uu — > tt, and eu — > et as well 
as the corresponding ones with the up quark replaced by a charm quark. We have calculated 
NLO and NNLO soft-gluon corrections to these cross sections. 

The Tevatron and HERA colliders provide a unique opportunity for precise study of FCNC 
couplings involving top-quark production. For some of these processes the uncertainty of the 
Born-level cross sections due to the choice of the QCD scale could be of the order of 100%. 
The search for FCNC single top-quark production at HERA is especially intriguing since a very 
recent experimental study of this process not only establishes a limit on the anomalous FCNC 
coupling but also presents an excess of events in the leptonic channel over the SM background 
(5 events observed while 1.31 ± 0.22 events are expected for SM background !) |2*lj . 

It is also worth mentioning that Tevatron and HERA are highly complementary in con- 
straining FCNC top-quark couplings: HERA has the best potential in limiting the k 7 coupling 
and is practically insensitive to Kz, while the Tevatron is sensitive to both kz and k 7 couplings. 
Therefore the combination of k 7 limit from HERA and Kz limit from the Tevatron will be the 
most stringent limit in k 1 - Kz plane. 

The main outcome of this paper is that we have found that higher-order QCD corrections 
stabilize the factorization and renormalization scale dependence of the FCNC cross sections 
down to the level of about 10% or less, and provide modest to significant enhancements to the 
leading order result. 
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